Particle number density fluctuations and pressure effects on structure formation 
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We provide a covariant and gauge- invariant approach to the question of how a first order pressure 
can be incorporated self-consistently in a cosmological scenario. The approximation is relevant, in 
the linear regime, to weakly self-interacting or warm dark matter models. We also derive number 
density fluctuations in which new modes appear because the number density fluctuations are no 
longer simply proportional to the density fluctuations. 
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I. INTRODUCTION 

<D ■ 

^ • According to current cosmological models, the large-scale structure, as seen today, formed from the evolution 
of gravitational instabilities which originated from inflationary processes in the primordial fluid. In studying the 
evolution of such instabilities, analytically in the fluid approximation^, one usually assumes a dust or radiation energy 
momentum tensor. The object of this paper is to consider some implications of modifying these assumptions. 
, First we note that while the two extreme fluid states are reasonable for most of the history of the universe, there 
■ are short, but at the same time important, specific periods in which matter is in the non-relativistic regime with a 

ON \ tiny although non-negligible pressure. These are periods during which collisions, or small random motions (velocity 
dispersion) of collisionless matter, give rise to a small but non-zero matter pressure. Even though in these cases the 
matter pressure in the background may be negligible, small disturbances (perturbations) of the spacetime induce a 
tiny pressure term which can modify the evolution of the inhomogeneities. Such a scenario may also have applications 
when studying weakly self-interacting matter (for details see Spergel and Steinhardt Qj) or warm dark matter (for 
more details of this approximation see Hogan and Dalcanton in the linear regime. These dark matter models have 
been the subject of much recent work using numerical integrations in the nonlinear regime, where they have been 
employed to investigate ways of resolving problems with the standard CDM model on small scales (see 111 ) in galaxy 
I ■ formation. 

Second we take the particulate nature of matter seriously. It must be born in mind that in structure formation one 
+3 . really needs to look at the clumping of a set of particles rather than at the evolution of energy density instabilities. 
Although in taking the dust assumption to be valid such a distinction is not relevant, this is not true for those 
periods in which matter is undergoing a transition to the non-relativistic regime and, consequently, for which the dust 
approximation is not accurate enough, even in the linear regime. To model the effects of the particulate nature of the 
matter we make use of thermodynamic arguments. An alternative Kinetic Theory approach is given in 

The method we use is based on the Ellis-Brum H covariant and gauge-invariant perturbation formalism, which 
has proved useful in different cosmological problems @, 0|, |§- It is one of the gauge-independent descriptions of 
relativistic perturbations, described in detail in the literature, see for instance ||, Q, ||. 

In our treatment of the fluid models we generalise the usual equations for the growth of cosmological perturbations 
for matter in a Friedman-Robertson- Walker (FRW) background to include first order pressure effects due to a non- 
relativistic component, in the radiation and in the matter dominated eras. Thus the simplifying hypothesis that the 
matter is dust, i.e. has vanishing pressure, is relaxed, in both cases, to include a small pressure due to collisions or 
thermal motions of the particles which make up the fluid. 
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x For high accuracy the full Boltzmann equations are used and these require numerical codes, which provide different insights 
into the process. 
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II. MATTER VERSUS ENERGY DENSITY PERTURBATIONS 



The starting point for the the Ellis-Bruni || formalism is the choice of a time-like congruence u a along which 
typical (comoving) observers move. Once this four-velocity has been selected all the physical quantities are split into 
their spatial and time-like parts using the spatial projector tensor h ab . The relevant physical quantities are defined 
as follows 

d a = a , 9 a = aD a 0, p a = a , v a = a , (1) 

p p n 

where D a = /i„Vf, is the spatial projection of the covariant derivative, 9 = u a a the expansion, p the pressure, p the 
energy density and a is a representative or average length scale defined by 3d/a = 9. In a FRW metric a is the 
background scale factor and we will use the notation 3H = 9 in the background. The quantities 5 a , p a , 9 a and v a 
describe the spatial inhomogeneity of quantities which are homogeneous and isotropic in the background. Each of 
them is spatial (contraction with u a vanishes), covariant and gauge- invariant to first order, as all of them vanish in 
the FRW background |^|. The scalar parts of perturbations described in ([!]) are obtained by taking the comoving 
spatial divergence, e.g., 8 = aD a 5 a - 

The general procedure to obtain evolution and constraint equations for the quantities defined in (|l|) starts from a 
suitable spatial projection of the Ricci and Bianchi identities as given in terms of the fluid variables once the field 
equations have been taken into account (see for example || for a pedagogical account). Here we shall follow the 
notation of Q. 

Restricting ourselves to a spatially flat FRW background, the equations governing the evolution of the comoving 
fractional density gradient 8 a and the comoving expansion gradient 9 a are (see |Hj equations 25,26) 



S a = 3wHS a -(l + w)9 a , (2) 
9 

— j 
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% - ^-H 2 p a + 2H9 a - aD a D b u b + ^L(p a + 3 Pa ) = 0, (3) 



where 



w = -. (4) 



Note that equation (|2|) omits second order terms proportional to a a bS b and Lo a b8 h where a a b and uj a b are the shear 
and vorticity tensors. 

We will now discuss perturbations in the particle number density. This discussion is motivated by the idea that a 
full understanding of large scale structure formation in the framework of gravitational instabilities requires an analysis 
of how particles, which comprise the cosmic fluid, clump together due to their mutual attraction. Although dealing 
with the energy density is the usual procedure, we believe that it is desirable to consider the particle number density 
n because, ultimately, particles are the building blocks of all the cosmic structures we see today. For dust, it does 
not make any difference whether one chooses p or n because no random motions are present and the internal energy 
vanishes identically, yielding p = mn, i.e. 

8 a = v a . (5) 

The situation is quite different when internal energy associated with thermal motions is introduced. The behaviour 
of v a requires a separate investigation. In B it is shown that from the Gibbs equation 



Tds 

it follows that 



(9 (6 » 



e a = S a - (1 + w)u a , (7) 

where 

anTU a s 

e a = > (8) 



2 



is a dimensionless entropy gradient. Perturbations which are non-dissipative (s = 0) and without spatial variation of 
the entropy (e a — 0) are called isentropic in Q because they have the same entropy at all the points of the spacetime. 
In the weaker case where no spatial variation of the entropy is allowed (e a — 0) the number density perturbations are 
related algebraically to the energy density perturbations by 

(9) 



" 1 + w 

In fact, e a — is a strong assumption which does not hold in general. Its departure from zero affects the growth of 
perturbations. From the particle conservation equation 



we get an evolution equation for v a 



n0 = O, (10) 



-Pa = 0. (11) 



1 



If we use the standard, but generally unphysical, equation of state w = constant then together with (g), the integral 
of (0) yields 

v a = T ~?— + K, k a = 0. (12) 



Thus, in general, when this simple equation of state is used, a new stationary (entropy) mode appears in the evolution 
of v a which could have an impact on the structure formation at certain stages of the evolution of inhomogeneities. 
The problem with (|l2| ) is that it makes sense only for the very special, but widely used, cases of dust and radiation. 
In more general circumstances such as in self interacting dark matter or where there are collisions (l) , the equations 
of state should involve two independent thcrmodynamical quantities, i.e. n and p. To deal with the inhomogeneities 
that arise in non-relativistic matter with a non- vanishing pressure we may follow two different approaches. 

(i) If we assume that the fluid has reached a collision-dominated equilibrium^ , we learn from kinetic theory that 
the equation of state in the non-relativistic regime is Jl0[ 

3 

p = mn+-p, (13) 

which together with the energy balance equation 

p + 6{p + p) = 0, (14) 
and the particle number balance equation (|l^) leads to 

P=-\o P , (15) 



with solution 



Also from ( |l3| ) we get the equation 



«(t) 5 - <»> 



a = Va + -zPa (!') 

p 2 



which, together with (Jill) and the conservation equations, determines v a . 



2 For an example of such a model see the paper by Spergel and Steinhardt [jij] where it is pointed out that under certain 
circumstances the cold dark matter behaves as a collision dominated gas 
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(ii) Without imposing the restriction of a collision-dominated equilibrium fluid, we could tackle the problem of 
a collisionless gas in which a small pressure arises because of velocity dispersion (random motions) of the particles 
making up the system^. In J3 it was shown that by neglecting higher powers of the velocity dispersion in the 
Boltzmann equation, the evolution equation for pressure is given by (|15|). We remark that in obtaining this equation 
no assumption about the equation of state was made. Only the absence of collisions and the smallness of the velocity 
dispersion were used. Despite the evident mathematical similarities to case (i), a crucial difference emerges here: for 
a collisionless gas the Gibbs equation (^|) does not apply and, as we shall show below, this leads to some new results. 

In the next section we discuss cosmological settings where the small pressure terms introduced above are included. 

III. NON-RELATIVISTIC MATTER PERTURBATIONS WITH PRESSURE 

In this section we investigate the role played by a first order pressure in the non-relativistic regime on an FRW 
background filled, in the first case, by a dust gas and, in the second, by dust and radiation. The point is that while 
the idealized background can be described by dust or dust and radiation, as soon as the dust is disturbed (by velocity 
dispersion or collisions) a small nonzero pressure arises. This means that the pressure is a first order quantity, as 
it vanishes in the background. This differs from the usual approach which takes an identically zero pressure in the 
background and in the real (perturbed) spacetime. We relax this assumption by allowing first order corrections to 
pressure due either to the fluid being in a collision-dominated regime (case (i) above) or through velocity dispersion 
without collisions (case (ii) above). This enables us to initiate a theoretical discussion of the effects being studied 
mainly qualitatively or numerically in the papers One-component and two-component models will be considered 
in turn. We note here that the perturbations will involve perturbations of the entropy. The precise relation of the 
covariant description of entropy perturbations to the metric description [ p"3| is not straightforward and will not be 
discussed here. 

A. Dust background 

For simplicity we use an Einstein-de Sitter background. As stated above we focus on perturbations for which 
pressure is a first order quantity. From the previous section we have that, both in case (i) and case (ii), p is given by 
the evolution equation ( |l5| ) leading to ( |l6| ) 



P = Po 



As p is now a first order quantity we get from (0) and 



(?)'■ m 



5 a - 2HS a + l^Sa + -V a D z p = 0, (19) 
2 p 



where the momentum balance equation 



{p + p)u a = -D a p, (20) 



has been used. Note that for large-scale perturbations the Laplacian term is negligible and we recover the standard 
evolution equation for dust. Hence, as expected, pressure plays a role only at small scales. 

B. Collision dominated fluid 

At this point we distinguish between the two cases discussed in section two. If we assume that the fluid is in 
collision-dominated equilibrium, then the Gibbs equation holds which, together with the equation of state (fill), leads 
to 

D a p-0, (21) 



3 Such a dark matter model has recently been described by Hogan and Dalcanton |jj 
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at first order, reducing (|T^) to the standard equation valid for dust at any scale. Expression ( pl| ) implies that the 
spatial gradient of p is a gauge-invariant second order quantity since it vanishes at first order. This allows us to go 
beyond first order and derive a second order equation for the pressure as we will now show. A fully covariant and 
gauge-invariant treatment for second order perturbations is not possible at the moment. The point is that in order to 
construct n-order gauge-invariant variables following Ellis-Bruni, we have to ensure that they are constant (usually 
zero) at all orders below n . In particular it is not clear how to find (if at all possible) a second-order meaningful 
variable for the perturbed density in such a way that it vanishes at zero and first orders. Bruni and co-workers provide 
a systematic method to tackle relativistic perturbations beyond the linear order, but one has to pay the price of losing 



the gauge-invariant character 12 



It is possible, however to obtain a second order gauge invariant constraint on the perturbations from the above 
results. Operate on equation ( |T5| ) with aD a /p and use the identity (p6|)f| below to obtain 



and then use, 



Pa+Pa- + ^0a- + % a = (22) 
p A p 6 



S a = -9 a + 0(2) 



p=-3\- P + 0(l) 



derived from equations (0) and (|14|) and where 0(n) indicates quantities of order n. Note that p a is of order 0(2). 



After a change in the independent variable the equations (22) become 



where a prime denotes a derivative with respect to a/ oq. According to the argument above and the fact that p and 
S a are first order quantities, equation ( p3| ) contains only second order quantities. It is one of the set of second order 
equations. We treat it as a constraint on the gauge invariant pressure perturbations in the sense that a full solution 
to all the second order perturbation equations will have to satisfy ( p^) although the actual solution may be more 
restricted. With this in mind we proceed to solve it on a dust background. In equation (|23] ) we use equations jl^ ) 
and d29|), below, for the dust background and the standard first order solution for dust 



5 a = Kt(-)+K-(^) V \ (24) 
\a J \a J 

where and K~ are constants of order 0(1), to obtain 

Here At, A2, and K~ are constants of order 0(1), so the products are of order 0(2) and K* are constants of 
order 0(2). The solution (^5|) is a simple and exact solution giving a limitation on the evolution of a second order 
quantity in the Ellis-Bruni formalism. The point is that we can apply the formalism because T) a p vanishes at the zero 
and at the first order approximation. 

C. Collision free fluid 

The more interesting and physically appealing situation is that of a collisionless gas for which the foregoing argu- 
ments are not applicable. Applying the operator aD a D 2 to ( |l5| ) and using the first order identity 

(aDafY - aDj, (26) 



Since p a is second order, the first order identity gives rise to second order terms and any corrections are of higher order. 
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where / denotes any first order scalar or tensor quantity, we get 



(aD Q D 2 p)' = -7 (^) aD a B 2 p, (27) 



a 



and hence 



aD a D 2 p = A° a (^y , A° a ^aB a B 2 p\ a=ao , A° a = 0. (28) 

For a dust background 



so equation (n9h reduces to 



which can be solved to give 



«. = /C|^] + /c(^ 3/2 -2#(^ (3 i) 



where K± = 0. Thus a new non-adiabatic mode is obtained which decays, albeit at a slower rate than the standard 
decaying mode for dust. Although the effect of the mode will be at small scales, as desired, the fact that it is decaying 
means that it is unlikely to be significant. 



D. Dust-radiation background 

The novelty in this subsection is that the dynamics of the background is given by a decoupled mixture of dust 
and radiation. We will consider perturbations of the matter component which again will be taken as dust in the 
background but with a non-vanishing first order pressure in the real spacetime. This situation mimics that of the 
cosmological fluid around the time of the transition from the radiation to the matter dominated era with a matter 
component which is acquiring a major role in the dynamics. In the following p r , p r and p m denote the zero order 
pressure and energy density of radiation and matter respectively, whereas p m is the first order pressure of matter as 
given by ©. 

Assuming that the fluids are decoupled and share the same 4- velocity^ we get 

{fir + Pr)u a + ^>aPr = 0, (32) 

p m u a + D a p m = 0. (33) 



From these equations and the equation of state for radiation we find that 

51 = Ap™, (34) 



Pr + Pr Pm 

where 



S a = a , p a =a . (35) 

Pr Pm 



5 Using the same velocity is valid when the two fluids are coupled. Here it is a simplifying assumption which enables us to 
provide an analytical treatment of a problem, that has only been treated numerically up to now, and we are able to obtain new 
results. The role of the assumption is made explicit, in particular in the comment following equation (pflj), where it is relevant. 
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Note that if we were to keep the matter component as dust both in the background and in the real spacetime, we 
would have ended up with u a = = D a Pr, i.e. radiation perturbations would not have been allowed. 

In order to get an evolution equation for 9 a we start with the Raychaudhuri equation 

6+^9 2 --D a u a + ^(p + 3p)=0, (36) 

where p = p r +p m is the total energy density and p — p r +p m is the total pressure. We have to be very careful in dealing 
with p since in the background p = p r but when we apply the spatial derivative operator we get D a p = D a p r + T> a p m . 
Bearing all this in mind the evolution equation for 9 a becomes 

6 a + 2H9 a = aB a B b u b - ^H 2 au a (l + w) - ^(Sp rP " a l + 3p m p™ + p m 5™), (37) 

where equation Q) has been used. 

Before dealing with the acceleration term we recall the identity 

D 2 (D a /) = D a (D 2 /) +\{P~ 3H 2 )D a f + 2/curl u a . (38) 
From this identity and bearing in mind that p m is a first order quantity, we have 



,-,/, . D a D p m D T> a p m , ... . 

aD a U ui, = —a = —a = — D p a . (39) 



It is important to note that the result obtained from commuting the operators D a and D 2 would have been different 
if we had used p r instead of p m . In the latter case a first order term ~ p r curl u a would have arisen. This means that 
consistency requires 

curl u a = 0, (40) 

i.e., in a radiation-matter decoupled mixture for which (a) both fluids share the same ^-velocity and (b) the matter 
pressure arises from perturbations so that p m is a first order quantity, curl w a vanishes at first order. 

Changing the independent variables from t to (a/ao), we get the evolution equation for 9 a , 



The background equations for dust plus radiation are 

3 /, 

— (1 + aa , j. ; , 

par pa 4 pa 
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p = Pm + pr = + P = Pr = ~r^, H 2 = —(1 + aa), (42) 



where 



Po a o ' Po a o ' 



From the evolution equation for p for a collisionless gas 
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(43) 



(44) 



and using the equation (]2q) we have 



which gives 



f: = -2Hp™ ^ (p™)' = --p-, (45) 



pZ = K n (^)\ Ka=0. (46) 
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Using similar reasoning for D 2 p™ and with the help of the identity 

(D 2 /)' = D 2 / — 2HD 2 f + /D a u a , 

from Q] , we get 

We now define the scalar parts of the perturbations by 

S = aD a S a , Z = aD a 6 a . 

From (W6T) it follows that 



M™ = 0. 



«D>'» = K 



K = 0, 



where if = aD a XI", and that 



where M = aD a M r a n . From equations (B3J) and (M) it follows that 



aD a [aD a (D b u b )} = -M (^)' 



M = 0, 



aD a (au a ) = -x(^)' 



Operating on equation (B7]) with aD a gives 



Z + 2#Z + 



l pm5 = _ M (^) 4 _ K (^) V + ^ m - |(1 + ^)i? 2 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 



(53) 



For the matter scalar perturbations, w = and equation (g) implies Z = —8. The evolution equation for 8 written in 
terms of the scaled independent variable 



x = aa 



Po a 



is 



" X + 2 -8'-l-±-,8 = F( x} , 



2x(x + l) 2x(x + l) 
where the prime now denotes a derivative with respect to the new independent variable x and 



F(x) = 



1+x 



M , 5 K fPoY 1 



Po x 2 \ Po J x 



is the inhomogeneous part. It is immediately apparent that 

8x = 3x + 2 



(54) 



(55) 



(56) 



is a solution to the homogeneous part of equation (|5J). The method of reduction of order then leads to a complete 
solution. The second solution of the homogeneous part is 



82 = 81 



dx 



{3x + 2) 2 x^/x~TT , 
which can be integrated using partial fractions to obtain 

3 



(3x4-2), 

5 2 = 2^ x + 1 + 4 ln 



VxTT- 1 



. -v/F+T + 1 



(57) 



(58) 
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The homogeneous form of equation (|54|) was obtained by Meszaros [|5| in 1974 who obtained a solution by transforming 
the equation to a hypergeometric form. Groth and Peebles jl6| obtained the solutions (|5^ , ^) and gave the asymptotic 
behaviour 

6 2 « 5Mr/8), x < 1, 
<5 2 « -8/(45r), ar > 1, 



where 



m\2 



v/3 (P'') 3/2 

The results are further described in Peebles's book [ pT| . 

The full solution to the inhomogeneous equation is given by 



1 

ESI 



S = C1S1 + C 2 S 2 + Si I —2 ( / E5\Fdx ) tfcc, (59) 



where _E = exp (/ (3a; + 2)[2x(x + l)]" 1 ^) = x^/x~TT. 
The expression ESiF can be written in the form 

^=7m{* + T + £ + £} («» 

where Ao = 9M/pq, Ai = ^j-K (p^/p^) . From this the particular integral can be calculated term by term. After 
integration, the first term i5 p i of the particular integral is given by 



Opl 



6i JikUwri dx ) dx 

= A + ^-(3x + 2) In ( ' 



2 v ' V3a; + 2 

The remaining terms are more difficult to integrate. From the first term it appears that the full inhomogeneous 
solution will contribute new modes. Among these it will contribute a constant mode and may contribute a mode 
similar to the second term in the solution 82 and hence modify the coefficient of the factor (3x + 2) which could affect 
the asymptotic behaviour. The significance of this has to be viewed with caution because the remaining, still to be 
determined terms, may contribute further modes or cause cancellations. From the physics one is led to expect that 
the contribution of a small pressure perturbation will only be decaying modes. 

Two other features of equation ( |54| ) are worthy of note. First, for large x the asymptotic form of the homogeneous 
equation is the usual equation for first order perturbations in dust. Second, even for initial conditions in which <5 and 
its first derivative with respect to x are zero, perturbations will arise from the influence of the function F(x). In other 
words small pressure effects or number density fluctuations can source density perturbations. This is possibly more 
significant, for the formation of structure, than the modification of the modes. 



IV. CONCLUSION 



For analytical as opposed to numerical modelling of the evolution of inhomogeneities in the universe, it is conven- 
tional to use a fluid approximation. At a detailed level this is at variance with reality because the matter is generally 
more particulate than hydrodynamics allows. Also most of the analytical literature assumes that the cosmological 
non-relativistic matter is dust. If we take the particulate nature into account, then at early times weak self-interactions 
or small random motions of collisionless matter give rise to a small but non-zero pressure in a fluid model. With this 
motivation we have discussed some of the implications of including a first order pressure in the matter distribution 
for fluid models. 

We have derived the relation between the energy density perturbations and the number density perturbations for 
the isentropic case and we show that in the simple non-isentropic case with p/ p — constant a new stationary mode 
appears. This new mode will affect the formation of structure at certain stages in the evolution of the inhomogeneity. 
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In the main body of the paper we derive an exact solution to the energy density perturbation equation when the 
matter pressure is non-zero at first order. The implications of assuming a non-zero first order matter pressure are 
quite deep as is illustrated in equation (p5|), which is an equation for second order gauge invariant quantities, and in 
equation ((3^), which together with ([32] ) shows that if D a Pm were equal to 0, then u a = = D a p r - This would reduce 
the problem to the standard case. 

The complete solution (|59|) is different from the usual solutions used in CDM approximations and from the Meszaros 
solution. Having such a solution may be useful as an analytical tool for understanding some features of recent nu- 
merical work on weakly self-interacting or warm dark matter models in the non-linear regime (![]. This is a subject of 
further investigation. 
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